Critical Behavior of Anisotropic Heisenberg Mixed-Spin Chains in a Field 
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We numerically investigate the critical behavior of the spin-(l, i) Heisenberg ferrimagnet with 
anisotropic exchange coupling in a magnetic field. A quantized magnetization plateau as a func- 
tion of the field, appearing at a third of the saturated magnetization, is stable over whole the 
antiferromagnetic coupling region. The plateau vanishes in the ferromagnetic coupling region via 
the Kosterlitz-Thouless transition. Comparing the quantum and classical magnetization curves, we 
elucidate what are essential quantum effects. 

PACS numbers: 75.10.Jm, 75.40Mg, 75.50.Gg, 75.40.Cx 



Quantized plateaux in magnetization curves as func- 
tions of a magnetic field for spin chains have been attract- 
ing much current interest. The trimerized spin-^ chain 



I, while the 
Li |, where m is 



exhibits a massive phase at m/m sa t = 
dimerized spin-1 chain at m/m sat = | 
the magnetization per unit period and m sa t is its satu- 
ration value. The presence of finite gap and plateau has 
further been discussed |^-^| and actually been observed 
fulfill for various polymerized spin chains and ladders. 
It may be the Lieb-Schultz-Mattis theorem jL2j and its 
generalization [^| Ji~^Jl^] in recent years that motivate such 
vigorous arguments. Oshikawa, Yamanaka, and Affleck 
(OYA) pointed out that quantized plateaux in magneti- 
zation curves may appear satisfying the condition 



integer , 



(1) 



where S is the sum of spins over all sites in the unit 
period. 

The OYA argument stimulates us to study quantum 
mixed-spin chains as well. An arbitrary alignment of 
alternating spins S and s in a magnetic field, which is 
described by the Hamiltonian 

N 

H = J2 [fa ■ «i)a + (Sj ■ S j+ l)a - H(S* + 8$\ , (2) 



with (S ■ s) a = S x s x + S v s y + aS z s z , shows ferrimag- 



netism |15 , instead of antiferromagnetism, and is an- 
other current topic from both theoretical p^-|2^] and ex- 
perimental psjffljfl points of view. As H increases from 
zero to the saturation field 



(S,s) — (1, |). Then a plateau may appear at m = i. 
At the Heisenberg point, the ground state of the Hamil- 
tonian (||) without field is a multiplct of spin (5 — s)N 
|l5| and thus has elementary excitations of two distinct 
types ]T^ , [l9| , ^l| . The ferromagnetic excitations, reduc- 
ing the ground-state magnetization, exhibit a gapless dis- 
persion relation, whereas the antiferromagnetic ones, en- 
hancing the ground-state magnetization, are gapped from 
the ground state. Therefore, at the isotropic point, m as 
a function of H should jump up to h just as the field is 
applied and remain unchanged until the field reaches the 
antiferromagnetic excitation gap 1.759 [ pl| . 
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FIG. 1. The ground-state magnetization curves for the 
quantum Hamiltonian (2) at a > 0, where QMC calculations 
at TV = 32 and T — 0.08 are also shown for comparison. 



H sat = a{S + s) + y/a 2 (S-s) 2 +ASs , 



(3) 



the OYA criterion (|l|) allows us to expect quantized 
plateaux at m = S + s — 1, S + s — 2, • • • 1 (or |). 
Since the low-energy physics of the model (g) is qual- 
itatively the same |l7],[l8|,^7j regardless of S and s as 
long as S ^ s, here, let us consider the simplest case 



Once we turn on the exchange-coupling anisotropy, the 
plateau is not so trivial any more. We show in Fig. [I] the 
zero-temperature magnetization curves of the anisotropic 
chains, which have been calculated by the numerical diag- 
onalization technique combined with a finite-size scaling 
analysis pQ| . In order to verify the reliability of our scal- 



1 



ing analysis, which is briefly explained later, we have car- 
ried out quantum Monte Carlo (QMC) calculations |3l]] 
as well at TV = 32 and N — 16, where, due to the small 
correlation length jl8],[l9| of the system, the data show no 
size dependence beyond the numerical uncertainty. Al- 
though the QMC findings are obtained at a sufficiently 
low but finite temperature, they fully suggest that the 
present diagonalization-based calculations well describe 
the thermodynamic-limit properties. As the model ap- 
proaches the Ising limit (a — ► oo), the plateau monoton- 
ically grows and ends up with a stepwise magnetization 
curve. On the other hand, the introduction of the XY- 
like coupling anisotropy reduces the plateau. Thus we 
take great interest in where and how the plateau vanishes. 
Alcaraz and Malvezzi JljJ showed that the ground state 
of the model ^ without field is in the critical phase over 
the whole region — 1 < a < 1. At the Heisenberg point, 
the model still lies in the massless phase but the low- 
energy dispersion as a function of momentum is quadratic 
pjj . For a > 1, the model is in the massive phase 
and its low-energy structure is well understood by the 
spin- wave dispersions 



,T - 



= \h 2 (s 



s) 2 -4S*scos 2 |to(S 



(4) 



which describe the sector of the magnetization J2j($j + 
s|) = M <{S- s)N and that of M > (S - s)N, respec- 
tively. Thus the introduction of the anisotropy essentially 
changes the nature of the model (||) and a fascinating 
physics must lie especially in the XK-like coupling re- 
gion. In the present article, we clarify how the quantized 
plateau of the (1, |) model behaves as a function of the 
anisotropy and aim to reveal critical phenomena inherent 
in quantum ferrimagnets. 

In order to investigate the quantum critical behavior, 
we carry out a scaling analysis on the numerically calcu- 
lated energy spectra of finite clusters up to N — 12. Let 
E(N, M) denote the lowest energy in the subspace with 
a fixed magnetization M for the Hamiltonian (||) with- 
out the Zeeman term. The upper and lower bounds of 
the field which induces the ground-state magnetization 
M are, respectively, given by 



+ (N,M) 



H 

H- (N, M) 



E(N, M -t 
E(N, M) 



1) — E(N, M), 
- E(N,M - 1). 



(5) 
(6) 



If the system is massive at the sector labeled M, 
H± (N, M) should exhibit exponential size corrections 
and result in different thermodynamic-limit values 
H±(m), which can precisely be estimated through the 
Shanks' extrapolation ]32|j . For the critical system, on 
the other hand, H± (N, M) are expected to converge to 
the same value as [B3f 



H±(N,M) ~ H(m) ± 



(7) 



where v s is the sound velocity and r\ is the critical index 
defined as {(Jq(j~) ~ (— l) r r~ n for the relevant spin op- 



erator a, which may here be an effective combination of 
S and s. 

Figure ^ was thus obtained, where we smoothly inter- 
polated the raw data H(m). The system is trivially gap- 
less at all the sectors of M for a < — 1 and should there- 
fore encounter a massive-massles phase transition in the 
AY"-like coupling region. Now we present in Fig. || the 
magnetization curves at a < so as to detect the transi- 
tion. Surprisingly, the plateau still exists at the XY point 
(a = 0) and the transition occurs in the ferromagnetic- 
coupling region. At an naive idea of relating the massive 
state with the staggered Neel-like order in the direction 
of the external field, we are never able to understand why 
the plateau is so stable against the AF-like anisotropy. 




FIG. 2. The same as Fig. 1 but at a < 0. 
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FIG. 3. (a) Scaled quantity N An versus a. (b) The central 
charge c and the critical exponent r\ versus a in the vicinity 
of the phase boundary. 
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FIG. 4. The linear extrapolation of the size-dependent 
fixed point a c (N, N + 2). 

The plateau length A N = H + {N, M) - H_(N, M) is a 
relevant order parameter to detect the phase boundary. 
When the system is critical, An should be proportional 
to 1/N because of the scaling relation (m). We plot in 
Fig. |^(a) the scaled quantity NAn as a function of a. 
NAn is almost independent of N in a finite range of a, 
rather than at a point, which exhibits an aspect of the 
Kosterlitz-Thouless(KT) transition Q. It is likely that 
the XF-like anisotropy induces the KT transition |35| ] 
followed by the gapless spin-fluid phase whose spin cor- 
relation shows a power-law decay. Let us evaluate the 
central charge c of the critical phase, which is expected 
to be unity. The asymptotic form of the ground-state 
energy pi 



E(N, M) 
N 



e{m) 



ircvg 



(8) 



allows us to extract c from the finite-cluster energy spec- 
trum provided v s is given. Here we calculate v s as 



N 
2^ 



[E kl (N,M)-E(N,M)} 



(9) 



where k\ = 2tt/N and Ek(N, M) is the lowest energy in 
the subspace specified by the momentum k as well as by 
the magnetization M. The size correction for the formula 
(0) is of order 0(1/N 2 ), which is essentially negligible in 
the present system. In Fig. ||(b) we plot c versus a and 
find that c approaches unity as the system goes toward 
the critical region. We further investigate the critical ex- 
ponent rj so as to verify the KT universality and to specify 
the phase boundary. In the critical region the asymptotic 
formula An ~ 2 r Kv s rj/N enables us to estimate rj. Since 
the KT transition holds rj = j at the phase boundary, we 
can evaluate the transition point a c from rj as a function 
of a. Figure ||(b) claims that a c = —0.41 ± 0.01, where 
c = 1.00 ± 0.01. The phenomenological renormalization- 
group (PRG) technique |3(| is another numerical tool to 
determine the phase boundary. Taking An as the or- 



der parameter, we extract the size-dependent fixed point 
a c (N, N + 2) from the PRG equation 



{N + 2)A N+2 (a) = NA N (a) 



(10) 



In Fig. | we plot a c (N,N + 2) as a function of 1/(JV+1), 
which is linearly extrapolated to a c = —0.57. The PRG 
estimate is somewhat discrepant from the above-obtained 
phase boundary. Here we should be reminded of Nomura- 
Okamoto's enlightening analysis p7| on usage of the PRG 
method. The PRG equation applied to a gapful-gapful 
phase transition yields a reliable estimate of the critical 
point, whereas, for a transition of KT type, the PRG es- 
timate is quite likely to miss the exact solution due to the 
incidental logarithmic size correction, encroaching upon 
the KT-phase region. Considering the limited availabil- 
ity of the PRG analysis, we may recognize the present 
PRG solution as the lower boundary of a c . 




FIG. 5. The ground-state magnetization curves for the 
classical Hamiltonian (2). 




FIG. 6. The ground-state spin configurations as functions 
of the field for the classical Hamiltonian (2) . 
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namely, the relevant correlation length is not divergent. 
The only point (a c , H c ) possesses the KT character. The 
long-lived plateau against the AY-like anisotropy, which 
is contrastive to the corresponding classical behavior, de- 
serves special remark and further investigation. We ex- 
pect magnetic measurements on anisotropic systems [p8[ . 

It is a pleasure to thank H.-J. Mikeska and U. 
Schollwock for helpful discussions. This work was sup- 
ported by the Japanese Ministry of Education, Science, 
and Culture through Grant-in-Aid No. 09740286 and 
by the Okayama Foundation for Science and Technology. 
The numerical computation was done in part using the 
facility of the Supercomputer Center, Institute for Solid 
State Physics, University of Tokyo. 

FIG. 7. The a-versus-H phase diagram for the quantum 
Hamiltonian (2). 



In order to elucidate how far from intuitive the present 
observation is, we compare it with the classical behavior. 
Let us consider the classical version of the Hamiltonian 
(|J) , where Sj and Sj are classical vectors of magnitude 1 
and i, respectively. We show in Fig. || the classical mag- 
netization curves and learn both similarity and difference 
between the quantum and classical behaviors. In the 
ferromagnetic and Ising-likc antiferromagnetic exchange- 
coupling regions they are quite alike, which is convincing 
in that quantum effects are supposed to be less significant 
in both the regions. However, the quantum behavior is 
qualitatively different from the classical one in the XY- 
like coupling region. The classical state of M — N/2, 
which is stable enough to form a plateau in the Ising-likc 
coupling region, is no more massive at a < 0.943. The 
spin configuration as a function of the field, revealed in 
Fig. [6| is suggestive in understanding the prompt col- 
lapse of the classical plateau with the increase of the 
Xy-like anisotropy. In the classical case, the spin con- 
figuration in the massive state of M = N/2 is stuck to 
(Sj,Sj) = (1, — 5). In other words, the plateau can not 
appear unless the configuration (Sj,Sj) = (1, — |) is re- 
alized. This is not the case for the quantum system. 
The quantum spin configuration in the massive state of 
M = N/2 generally depends on a and exhibits a quan- 
tum reduction from the classical Neel-like state. At the 
Heisenberg point, for example, the quantum averages of 
the sublattice magnetizations per unit cell in the massive 
state are estimated to be 0.793 and —0.293, respectively. 
It must be the quantum spin reduction that makes the 
massive state tough against the AF-like anisotropy. 

In sum the quantum mixed-spin Heisenberg model (j^) 
with (S, s) — (1, i) shows the three distinct phases at the 
sector of M = N/2; the plateau phase, the gapless spin- 
fluid phase, and the ferromagnetically ordered phase, as 
illustrated in Fig. [7]. The plateau appears for a > —0.41, 
including a ferromagnetic-coupling region. We note that 
on the boundary of the plateau phase except for the point 
(a,H) = (a c ,H c ) = (0.41,0.293) which is indicated as 
KT point in Fig. m, the plateau length is generally finite, 
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